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Exchange  and  Correlation  Effects  in  an  Inhomogeneous  Electron  Gas* 

W.  Kohn  and  L.  J.  Sham 
University  of  California,  San  Diego 
La  Jolla,  California 


ABSTRACT 

From  a  theory  of  Hohenborg  and  Kohn,  approximation  methods 
for  treating  an  inhomogeneous  system  of  interacting  electrons  are 
developed.  They  lead  to  self-consistent  equations  analagous  to  the 
Hartree  and  Hartree-Fock  equations,  respectively.  In  these  equa¬ 
tions  the  exchange  and  correlation  portions  of  the  chemical  poten¬ 
tial  of  a  uniform  electron  gas  appear  as  corrections  to  the  electro¬ 
static  potential. 
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Exc  hange  and  Correlation  Effec is  in  an  Inhomogeneous  Electron  Gas* 


W.  Kohn  ind  L.  J.  Sham 
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La  Jolla,  California 


It  has  been  shown1  that  the  ground  state  energy  of  an  interact¬ 
ing  inhomogeneous  electron  gas  in  a  static  potential  v(r)  can  be  written 
in  tne  form 

n(r )  n(r'  ) 

E  =  f  v(r)  n(r)  dr  +  ^  - dp  dr*  +  G[n]  (l) 

J  ~  ■-  >jJ  I  r*  _  -r  1  i 


where  n(r)  is  the  density  and  Gjn]  is  a  universal  functional  of  the 
density.  This  expression,  furthermore,  is  a  minimum  for  the  correct 
density  function  n(r).  In  this  note  we  propose  an  approximation  for 
Gfn.,  which  leads  to  a  scheme  analogous  to  Hartree's  method  but  con¬ 
tains  the  major  part  of  the  effects  of  exchange  and  correlation. 

We  first  write 


G[n]  h  T  [n]  +  E  fn] 


xc‘ 


(2) 


where  T  [n]  is  the  kinetic  energy  of  a  system  of  non-interacting 
s 
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electrons  with  density  n(r)  and  E  [n]  is,  by  our  definition,  the 

xc 

exchtnge  and  correlation  energy  of  an  interacting  system  with  density 
n(r).  For  an  arbitrary  n(r)  one  can  of  course  give  no  simple  exact 
expression  for  E  [n].  However,  if  n(r)  is  sufficiently  slowly  vary- 
ing,  we  can  show^  that 

ExJn]  =  J  n(r)  exc(n(r))  dr 


(3) 
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where  e  (n)  is  the  exchange  and  correlation  energy  per  electron  of 

vn  '  ' 


xc 


a  uniform  electron  gas  of  density  n.  We  propose  (3)  as  our  approxi¬ 


mation  even  when  n(r)  is  not  slowly  varying  and  shall  regard  e  as 

XC 
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known  from  theories  of  the  homogeneous  electron  gas. 

From  the  stationary  property  of  (l)  we  now  obtain,  subject  to 
the  condition 

6n(r)ar  =  0  ,  (4) 

v 

the  equation 


: &n(r) 


here 


6T  rn] 


XC 


n 


(r) 


dr  =  0  ; 


(5) 


<?(r)  =  v(r) 


I  r-r  ’ 


dr 


(6) 


and 


U  (n)  =  ^—  (ne  (n) )  (7) 

xc  '  an  \  xc  / 

in  the  exchange  and  correlation  contribution  to  the  chemical  potential 
of  a  uniform  gas  of  density  n. 

Eqs.  (4)  and  (5)  are  precisely  the  same  as  one  obtains  from  the 
theory  of  ref.  1  when  applied  to  a  system  of  non-interacting  electrons, 
moving  in  the  given  potential  c?(r)  +  ^n(r)j  .  Therefore,  for  given 
cp  and  p,  one  obtains  the  n(r)  which  satisfies  these  equations  simply  by 
solving  the  one-particle  Sc'nroeainger  equation 


'I 


and  setting 


N 

n(r)  Pife)|2, 

i  =  1 


(9) 


'.mere  N  is  the  number  of  electrons. 

Eqs.  (6)  -  (9)  have  to  he  solved  self-consistently :  One  begins 
with  an  assumed  n(r),  construct.;  cp(r)  from  (6)  and  p  from  (j),  and 
finds  a  new  n(r)  from  (8)  and  (9).  The  energy  is  given  by 


e . 

l 


n(r)  n(r'  ) 


dr  dr '  + 

[  n(r ) 

[«  (n(r)\ 

1  -  p  fn(r )) 

J 

xc  \  v~V 

xc  \  ~  /J 

(10) 


II  ic  physically  very  o;.ti6factory  that  p.  appears  in  Eq.  (8) 

Av* 

us  an  additional  effective  potential  so  that  differences  of  p  lead 

xc 

to  pressures  on  the  electron  fluid  in  a  manner  familiar  from  thermo- 
ci.yna.nics . 

The  results  of  this  procedure  are  exact  in  the  following  limit¬ 
ing  cases:  (l)  when  n(r)  is  slowly  varying,  (2)  when  cp  dominates  over 
P  (as  near  an  atomic  nucleus),  and  (3)  when  n  is  very  large. 

XC 
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The  method  differs  from  Slater's  method  of  the  "exchange  hole" 

in  two  ways:  (l)  Slater's  method  does  not  include  correlation  effects,^ 

end  (2)  more  importantly,  the  Slater  effective  exchange  potential  is  3/1 

of  P_  ,  the  exchange  part  of  the  correction  potential  p  in  our  method. 

X  xc 

One  can  verify  by  a  study  of  the  slowly  varying  case  that  our  result  is 
correct. 


It  is  also  possible  to  obtain  a  scheme  which  includes  exchange 


effects  exactly,  by  writing  in  place  of  Eq.  (3) 


•  1 

ail 

x  ■  1 

n 

-  E 

X 

n 

v.hcr 


n 


*  l  J 

tens  ty  n(r) 

equa  -ion. 


is  the 
This 


+  y"dr  n  ( 

t  \ 

•c  (“■ 

(r) 

0 

exchange 

era 

tgy  in 

a 

Hartree-Fock  system  with 

leans  in 

pi. 

JO  of  1 

(8) 

1  to  a  Hartree-Fock  type 

(11) 


The  same  approach,  when 
resuats  for  a  number  of  other  el 
state  applications  are  cohesive 
laces,  electronic  specific  heat, 


appropriately  extended,  also  gives  simple 
-•e ironic  properties.  A* r.ong  the  solid 
.energies,  elastic  constants,  Fermi  sur- 
arid  energy  bands.  Details  will  be  given 


in  a  future  publication. 
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